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INTRODUCTION 

Sonoluminescence ("light from sound”) is the result 
of extremely nonlinear pulsations of gas/vapor bubbles 
in liquids when subject to sufficiently high amplitude 
acoustic pressures. In a single collapse, a bubble’s 
volume can be compressed more than a thousand- fold 
in the span of less than a microsecond. Even the sim- 
plest consideration of the thermodynamics yields pres- 
sures on the order of 10,000 ATM. and temperatures 
of at least 10,000K. On the face of things, it is not 
surprising that light should be emitted from such an 
extreme process. For a general review of some of the 
remarkable experimental results issuing from the study 
of sonoluminescence in acoustically levitated bubbles 
wc refer to a number of articles [Roy 1994; Crum 
1994; Crum & Roy 1994; Putterman 19951 The phe- 
nomenon has achieved prominence in the public eye 
[Browne, 1996; Glanz, 1996; Knight, 1996]. The in- 
crease in scientific interest can be gauged by noting 
that, while as recently as 1990 (the year after Gaitan 
discovered light from a single bubble) there were 4 
SL-rclated articles in refereed journals, this number 
increased to 13 in 1995, 24 in 1996, 31 in 1997, and 
14 refereed journal articles have already appeared in 
1998! 

Experiments [Gaitan et al. 1992a,b, 1997; Barber et 
al. 1991, 1992, 1994; Hiller et al. 1992, 1994; Holt el 
al. 1994, 1996; Lofstedt el al. 1995; Matula et al. 1995, 
1996, 1997a, b,c,d; Lepoint el al. 1997, among others] 
have uncovered the existence of many unexplained 
features surrounding the observation of sonolumines- 
cence. Of these, we have identified four which are 
fundamental, and in which gravity plays (or is pre- 
dicted by theory to play) a critical role. These are: 

• the light emission mechanism itself, 

• the mechanism for anomalous mass flux 
stability, 

• the disappearance of the bubble at some 
critical acoustic pressure, and 

• the appearance of quasiperiodic and cha- 
otic oscillations in the flash timing. 

Gravity, in the context of buoyancy, is implicated 
in all four of these unexplained phenomena. 

PROBLEM STATEMENT 

The scientific objectives of the planned research 
are: 


(1) To test the predictions of the fractolumines- 
cence light emission theory of single bubble sonolumi- 
nescence [Prosperetti 1997] which relies on the asym- 
metry induced by buoyancy-induced translatory oscil- 
lations to postulate that a jet forms during the final 
stages of bubble collapse. This hypothesized jet 
shoots through the bubble interior and makes contact 
with opposite bubble wall at near Mach speeds and 
'fractures’ the water, with resultant light emission. 
Vortices in the liquid carry off excess linear momen- 
tum [Longuet-Higgins 1996], so that the process can 
repeat the next acoustic cycle. This theory has many 
specific predictions, which can be summarized here by 
noting that, if gravity is reduced to the ambient orbital 
level of 10 5 g, then this theory would predict a com- 
pletely different picture of the parameter space and its 
regions of stability and light emission than the one 
which has recently been uncovered in 1 g experimenta- 
tion [Holt and Gaitan, 1996, 1996a; Gaitan and Holt, 
1997]. Simply by making careful measurements of the 
details of light emission and mechanical oscillations of 
bubbles in pg, where buoyancy-induced coupling of 
translation and oscillation motions is absent, we can 
definitively establish the viability of this theory. 

(2) To measure the mass flux and mechanical sta- 
bility boundaries which constrain the appearance of 
light emission in the parameter space of acoustic pres- 
sure, dissolved gas concentration and equilibrium ra- 
dius. These measurements will test both the fractolu- 
minescence and the chemical reaction/diffusion theo- 
ries that model the non-diffusivc anomalous mass flux 
observed in lg experiments [Lofstedt et al. 1995; Holt 
& Gaitan 1996a]. The latter theory [Lohse et al. 
1996a,b] relies on the assumption that gravitational 
effects are not important in the convective-diffusion 
treatment of mass transfer. 

(3) To measure the precise values of acoustic pres- 
sure and equilibrium radius where a light-emitting 
bubble disappears, an unexplained phenomenon which 
is a ubiquitous feature of lg experiments. The disap- 
pearance represents an upper bound on the acoustic 
pressure (and hence energy input to the system) which 
can be employed to drive the bubble oscillations. 
These measurements will directly test a theory [Cordry 
1996; Matula et al. 1997c] which postulates that the 
variation in buoyant force with each acoustic cycle's 
expansion and contraction of the bubble leads to a loss 
of levitation above some critical acoustic pressure. 

(4) To test whether chaotic and quasiperiodic tim- 
ing of the flashes observed in lg [Holt et al. 1994] is 
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due to buoyancy-induced effects or germane to the 
phenomenon of light emission itself. The nonlinear 
coupling of translational and radial motion via the 
buoyant force is postulated to result in quasiperiodic 
and chaotic modulation of the volume oscillations, but 
this would not occur in (ig. Simultaneous monitoring 
of the timing of individual flashes during detuning of 
the system will allow us to determine the presence or 
absence of the nonlinear effects, and thus to prove or 
disprove the conjecture. 

BACKGROUND 

Rather than giving a comprehensive review of light 
emission and nonlinear bubble dynamics in general, 
we will use the four above features to both motivate 
and organize our background section. Reviews of 
nonlinear bubble dynamics are numerous; a compila- 
tion is available in Leighton [1994]. Since the re- 
search is so topical, no comprehensive review of com- 
peting theories for the light emission mechanism can 
be written. We briefly list some of the top contenders, 
which include electron conduction-quenched plasma 
[Moss et ah, 1997], shock-wave mediated 
Bremmstrahlung [Wu & Roberts 1994], plasma dis- 
charge [Lepoint et al. 1997], superradiance [Trenta- 
lange & Pandey 1995], acoustic superresonance 
[Brenner et al. 1996], quantum vacuum fluctuations 
[Eberlein 1995], confinement of post-ionization elec- 
trons [Bernstein and Zakin 1995] and collision- 
induced dipole emission [Frommhold & Atchley 
1994]. Other than the fractoluminescencc theory dis- 
cussed above, only one of these various theories for 
the light emission mechanism would be directly influ- 
enced by gravity, and that is the plasma discharge the- 
ory, which relies on a shape instability (Rayleigh- 
Taylor mechanism) to concentrate electric charge on 
the surface of a critical microjet that penetrates the 
collapsing bubble interior. The net result is "lightning" 
inside the bubble due to plasma discharge associated 
with spray electrification of the jet as it breaks up into 
droplets. While this theory has not gained widespread 
acceptance, it nevertheless has been able to explain 
some features of the light emission. Except to say that 
it remains a contender and falls under the heading of 
"affected by gravity", we will not discuss this theory 
further. In the next section, we concentrate on one 
theory in particular which has the potential to explain 
in a unified fashion many of the experimental obser- 
vations. 

Time-varying buoyancy: the direct effect of gravity 
Consider the situation of a single air bubble in 
water, subject to a harmonic acoustic field which is 
also a standing wave. The acoustic wavelength is 
much larger than the bubble radius R 0 ( kR 0 is small, 
where k is the acoustic wavenumber). On the acoustic 
time scale (typical acoustic periods are 50 microsec- 
onds), the bubble will contract during the compression 
phase of the field. The pressure and temperature will 


increase. During the expansion phase of the external 
field, the bubble expands. Thus the bubble volume, 
V(t), can change significantly during an acoustic cycle. 

An oscillating bubble in a standing pressure wave 
experiences a nonlinear body force dominated by the 
large relative compressibility between the bubble and 
the host fluid. The force depends on the time-average 
of the product of the bubble volume and the pressure 
gradient, known as the primary Bjerknes force [Ray- 
leigh 1894; Bjerknes 1906; Yosioka & Kawasima 
1955; Doinikov 1994]. Its functional form is F a < 
V(t)*V P a (r,t) >, where <...> denotes a time-average 
over one acoustic cycle, and PJr,t) is the acoustic 
pressure field. F a acts in the direction of the pressure 
gradient of the standing wave, and thus it can balance 
the buoyant force and cause stable trapping of a bub- 
ble. 

Finally, the bubble experiences a buoyant force F b 
= p gV(t) directed opposite the acceleration due to 
gravity. In practice, the fact that this force varies in 
time with the bubble volume is usually ignored, due to 
the small change in volume coupled with the inertia 
and drag resistance to translational motion in the liq- 
uid. Thus, one normally considers an acoustically 
levitated bubble to have a fixed vertical position at 
constant acoustic pressure which will depend on the 
ratio of the magnitudes of the forces Fi/F a . However, 
in the parameter regimes where sonoluminescence 
occurs, the ratio of can reach 10 3 and greater, 

as shown in Fig. 1 . 



Figure 1 . The measured ratio of V ^ / V 0 for a bubble as 
a function of pressure. Light emission begins very near the 
peak value of the normalized response. The dissolved air 
concentration in the water was 15% of the saturation con- 
centration at 1 atmosphere ambient pressure. The equilib- 
rium or ambient radius varies with pressure in a complicated 
fashion, and is implicit in these measurements. Derived 
from Gaitan and Holt (1998). 

The resulting ratio of FJF a can vary by a factor O(10 3 ) 
during a single acoustic cycle — thus the bubble must 
undergo translatory oscillations driven by the time- 
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varying buoyant force F b \ This is a fundamental re- 
sult, independent of any proposed theory for sonolu- 
minescence. It is the basis for the fact that gravity 
plays a tremendously important role in these nonlinear 
bubble oscillations (despite that fact that a straightfor- 
ward calculation of a static Bond number for earth- 
based levitation yields p g(Ro? / 2 <7 ~ 0(10 7 ) for a 

10 micron bubble!). The effect of gravity lies in the 
magnitude of the changes in F b . The ratio F\/F a can be 
considered a coupling strength between the radius of 
the bubble and the position of the bubble’s center of 
mass. Because of this nonlinear coupling, volume os- 
cillations can drive translatory oscillations. If gravity 
is effectively removed, this coupling constant will be- 
come much smaller than unity, and thus its effects will 
virtually disappear. 

INVESTIGATIVE APPROACH 

Our short-term plan is to perform lg laboratory, 
drop tower, and parabolic flight experiments. Certain 
basic questions can be answered scientifically by drop 
tower and parabolic flight tests. But going beyond a 
"light or no light” level of sophistication will require 
the experimentation time to change key parameters. 
Significant experimentation, as we have documented, 
has been performed in lg leading to a rationale for an 
experimental design in which the critical measurables 
are already identified and the relevant parameters are 
well characterized. Since Og experimentation time is 
limited, the scope of our plan is finely focused: this 
investigation is mainly concerned with probing 
mechanisms for the light emission and related me- 
chanical stability, not with a parameter study. In the 
section that follows we detail the experimental plan; 
here, we present the plan’s essential features for rapid 
consideration. The crucial host fluid is water. The 
important dissolved gases are air and argon. Con- 
trolled dissolved gas concentration must be achieved in 
preparation and/or storage. The acoustic system must 
be capable of generating controlled acoustic pressures 
ranging from 1/10 ATM to 3 ATM. The dissolved gas 
concentration must be varied, either by pre-flight 
preparation, in-flight preparation, or variation of the 
test cell’s ambient pressure. Temperature must be 
controlled. The critical measurables and their associ- 
ated method are: 

• acoustic pressure versus time, P a ( t) 

(in situ hydrophone) 

• instantaneous radius versus time, R( t) 

(laser light scattering) 

• maximum and equilibrium radius, & Ro 

(imaging) 

• flash timing, intensity and duration 

(time-resolved photometry) 

• flash spatial distribution 

(multiple coincidence detectors) 


We emphasize that these are all techniques with 
which we have extensive experience from ground- 
based laboratory experimental work. This is not to say 
that the experiments have all been done, but that the 
measurements we propose are eminently feasible. 

Method and Diagnostics 

The experimental techniques will build on an 
optical scattering/acoustic levitation/imaging system 
used in Holt and Gaitan [1996, 1996a, 1997] and a 
photometric system used by Holt el al. [1994], and by 
Matula & Roy et al. [1997a,b]. This section will sim- 
ply highlight features of the apparatus and method, as 
well as discuss calibration and procedures. 

Figure 2. Schematic of the experimental setup. 



To monitor the oscillations of the bubble wall, it is 
illuminated by a linearly-polarized laser beam, and the 
scattered intensity as a function of time at a specific 
angle from the forward is recorded. The voltage out- 
put of a photodetector (typically a fast photodiode, but 
PMT’s are useful for low intensity work) positioned at 
the aforementioned angle is the time signal of interest. 
For purely spherical volume oscillations, a monotonic 
transfer function can be theoretically derived and ex- 
perimentally calibrated [Holt & Crum 1990] relating 
the voltage and the bubble radius R(t). 


To determine the presence and nature of shape os- 
cillations, a variety of techniques will be utilized in the 
manner of Holt & Gaitan [1996]. First, to detect the 
onset and frequency of shape oscillations, the scattered 
light intensity V exp (t) can be monitored for either a pe- 
riod-doubling or a rapid amplitude change [Gaitan & 
Holt, 1996]. Time-gating techniques can be used to 
provide a trigger criterion for shape oscillations cou- 
pled to the ringing oscillation after the main bubble 
maximum expansion. The values of acoustic pressure 
P a and equilibrium radius R 0 at the onset of shape os- 
cillations will be measured. Secondly, to determine 
the mode of oscillation, high-resolution, shorl- 
illumination-pulse video will be used. 


High-resolution video imaging with short (clpsec) 
pulse illumination will be employed to obtain me- 
chanical response information, Ro, Rmax, shape mode, 
and translational motion can all be obtained via this 
technique by improving on the system used by Holt 
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and Gaitan [1996]. Digital image processing tech- 
niques will be utilized to derive statistics from individ- 
ual images. 

Calibration 

1 . Cell pressure: 

The method used to calibrate the cell pressure 
is outlined in Crum [1983]. A small hydro-phone will 
be permanently mounted inside the cell. The pressure 
gradient in the gravitational direction will be meas- 
ured. While levitating a bubble in the cell in lg we 
measure its equilibrium position with respect to the 
measured pressure gradient in the cell, its equilibrium 
radius, and the hydrophone voltage as the pressure is 
varied. Using the equation expressing the balance of 
the acoustic force and the buoyant force for small os- 
cillations to solve for the pressure, a calibration con- 
stant can be determined from fitting techniques. It is a 
null method, and is very sensitive to small changes in 
pressure. As a check, numerical fitting of the equa- 
tions of motion for small -amplitude oscillations will be 
performed with the calibrated photodetector output as 
the data to be fit by varying the acoustic pressure am- 
plitude in the model. Gaitan & Holt [1998] have veri- 
fied that, while the Rayleigh-Plesset equation cannot 
predict where in the P& R space the asymptotic equi- 
libria exist, it can accurately predict one of the triple 
(P a , Rq, Rmax) given the other two. Thus we can cali- 
brate in situ in Og as well. 

2. Bubble radius 

The quantity I exp (R,t), the incident scattered 
light intensity falling on the photodetector, is equal to 
lo J J l re i (R(t)) d6 dtp , where the limits on the 
spherical angle variables 0 and (p in the integral are 
determined by the particular photodetection scheme 
used. I 0 is the intensity incident on the bubble. I re i is 
the normalized component of the intensity (with linear 
polarization parallel to the scattering plane) scattered 
into an angle 6 from the forward. I rei is calculable 
from Mic theory. Thus, the photodetector output 
V exp (R,t) - F / J l rel (R(t)) dG dtp , where F is an 
apparatus-dependent constant to be determined empiri- 
cally. 

First, a stably oscillating bubble is obtained in the 
cell and moved into the beam. The output voltage is 
monitored on an oscilloscope, and the time average is 
recorded. During calibration, only linearly oscillating 
bubbles are used, since R 0 = <R(t)> only for linear 
oscillations, and hence <I exp (R,t)> = hxp(Ro) only for 
linear oscillations. As soon as V exp (R 0 ) has been re- 
corded, an independent measurement of R 0 is made by 
instantaneously turning the sound field off while ac- 
quiring video images [Holt & Gaitan, 1996a]. The 
video image can be calibrated by a known scale. The 
calibration constant F is then determined by taking the 
ratio of the experimental voltage to the relative inten- 
sity Dividing the photodiode output by F gives the 


experimental relative intensity, which can then be used 
to find the radius. In this fashion, bubble radii can be 
obtained to within about 4% [Holt 1989]. 

An alternate technique has recently been intro- 
duced [Lentz et al. 1995]. This technique relies on the 
angular variation in the Mie scattering, where the 
number and location of intensity maxima are unique 
for a given bubble radius (at fixed laser wavelength). 
Thus, an absolute determination of R is possible. 
While the technique is very accurate, it is much slower 
than the previous technique. Nevertheless, it will be a 
valuable independent check on the fixed-angle method, 
and is able to provide in situ calibration without ac- 
cessing the interior of the acoustic test cell. 

Procedure 

We present here the generic procedures for 
obtaining data by varying the accessible para-meters of 
the system. These can be divided into two classes: 

1 , Mechanical parameters 

P a will be varied between 0.5 and 1.5 ATM, 
depending on the dissolved gas concentration, and 
increasing to higher pressures if the pressure for disap- 
pearance is not the same as in lg. The practical prob- 
lem this introduces in lg, i.e. changing the position of 
levitation, can be addressed with feedback control on 
the vertical translational stage via a PC-controlled 
stepper mike motor. Pressure control itself will be 
accomplished via a PC-controlled step sequence, in- 
creasing and decreasing the pressure in a controlled 
fashion with operator intervention possible at every 
step. 

R 0 will be implicitly varied via the constraints of 
mass flux equilibrium and shape instability. From our 
experience in lg, the approximate range should be 1 - 
25 microns in equilibrium radius, and 20 to 70 microns 
in maximum radius. This may, of course, change dra- 
matically with the removal of buoyancy-induced 
translatory oscillations. 

The acoustic frequency /^will be fixed at the reso- 
nance of the acoustic cell; with temperature control 
this should remain constant. When investigating 
variations in flash timing, /</ will be varied no more 
than a few hundred Hz from the resonance frequency. 

2. Material parameters 

The host liquid for all experiments will be 
water, distilled, micro-filtered and de ionized to 
1 7Mohm resistivity at 20°C. Two gases will be used, 
air and argon. They will be dissolved in water at 
varying percentages of saturation with respect to 100% 
saturation at 760 mm-Hg ambient pressure at 20°C for 
air over water. Air as a mixture of gases in the pro- 
portions found at sea level will be dissolved in equili- 
bration at ambient pressures from 400 mm Hg down to 
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the vapor pressure. Argon will be prepared in the 
same way, only the absolute pressure of the argon will 
reflect its partial pressure in air at sea level. Thus, 
argon will be dissolved in equilibration at ambient 
pressures from 8mm down to 3mm Hg. 

Data Acquisition 

The optical Mie scattering data will be the 
primary source for both information and triggering 
purposes for the other data acquisition methods. A 
photodiode or photomultiplier tube with a scattering 
laser-line bandpass filter will detect the scattered light, 
A calibrated system will yield both equilibrium and 
time-varying radius values. This signal will be digit- 
ized at 2 GS/s, and minimum 10 bit resolution, with 
record lengths up to 1Mb. 

The Mie signal will be used to trigger a short light 
pulse (diode laser, superbright LED or short halogen 
strobe light) backlighting the bubble for obtaining high 
resolution video/CCD images Appropriate optical 
filters will be used to block out the Mie signal. 

Matched photomultiplier tubes (or a mini-array) 
will be used to detect the sonoluminescencc flashes. 
Flash timing circuitry will include a Constant Fraction 
Discriminator and Time-Amplitude Converter. Rela- 
tive intensity measurements will be made using pulse 
height analysis with appropriate coincidence criteria, 
both to yield intensity as a function of parameter 
variation, and to detect spatial anisotropies in the 
emitted light. 

MICROGRAVITY RATIONALE 

The aspect of gravity which affects the proposed 
measurements is the buoyant force on the gas bubble. 
This is unavoidable in a material sense, since the den- 
sity and compressibility contrast are required to obtain 
volume oscillations in an isotropic pressure field. In 
addition to the transla-tional oscillations discussed in 
the previous sections, the presence of the buoyant 
force (and the need to overcome it) results in: 1) a 
perturbation from the ideal spherical symmetry in 
which one would like to perform experiments; 2) a 
requirement of a minimum pressure below which ex- 
peri-ments cannot be performed simply because the 
bubbles will rise to the container surface; 3) non-linear 
acoustic streaming, which will affect measurements of 
mass flux and mechanical stability. 

We can quantify our discussion in the background 
section by considering the difference in the nonlinear 
coupling constant Fi/F a in lg and pg. To obtain sim- 
ple order-of-magnitude results in this section we ig- 
nore the dependence on the relative phase of the bub- 
ble oscillation with respect to the acoustic field. The 
time-average ratio FJF a is exactly unity on earth when 
the bubble is displaced above the central antinode of 


the standing wave acoustic field to a point where the 
gradient of the field is strong enough so that the forces 
balance — this is the definition of acoustic levitation. 
As long as the volume oscillations of the bubble re- 
main on the order of 25% or so, then the volume in the 
force equations V(t) can be replaced by V 0 , the equilib- 
rium volume of the bubble. As is borne out in labora- 
tory results, for these mild oscillations the bubble re- 
mains at the location given by Ff/F a = 1 ~ p g f 
K(P a kfz , given suitable low -amplitude approxima- 
tions, where K is a material constant, P is the acoustic 
pressure amplitude, k is the acoustic wavenumber and 
z is the vertical height above the antinode. Thus z = 
p g / K(P a kf has no time dependence and depends 
only on P a for low amplitudes. 

However, as Fig. 1 shows, the ratio V max /V 0 f° r the 
high-amplitude oscillations during sonoluminescencc 
can exceed 10 3 . Thus in lg the ratio of Ff/F a can vary 
roughly by a factor / <V(t)> = 1 0 3 during a single 
acoustic cycle. This means that the gravity -induced 
nonlinear coupling to the translational mode will be a 
maximum precisely when sonoluminescencc occurs! 

If, however, gravity is reduced to orbital ambient 
levels of < 10' 4 g, this means that F/F a is also reduced 
to 0(1 O' 4 ) or less, since F b is linear in g. Thus even in 
the sonoluminescence regime, the coupling strength is 
only of order 10' 3 . A bubble in microgravity will thus 
be located almost precisely at the acoustic pressure 
antinode, where the field gradient is zero. The nonlin- 
ear coupling between the volume mode and translation 
mode is removed. 

The time scales in the problem dictate the need for 
extended periods of low acceleration. There are, of 
course, very short lime scales involved in the problem, 
such as the picosecond scale for individual sonolumi- 
nescence flash widths. However, beginning with the 
diffusive time scale for the growth, dissolution or 
equilibration of bubbles (on the order of seconds), the 
need for long time experiments become clear. For a 
given dissolved gas concentration, both mechanical 
and photometric data must be gathered as phenomena 
develop on the diffusive time scale at a constant 
acoustic pressure. Then the pressure must be incre- 
mented and measurements repeated while still at the 
same dissolved gas concentration, until the critical 
pressure and radius where the bubble disappears. Such 
a sequence, even if fully automated, could require tens 
of minutes. 
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